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We present a numerical study of the space-space and time-time components of the gluon propagator at equal 
time in the minimal Coulomb gauge, and of the gluon and ghost propagators in the minimal Landau gauge. This 
work allows a non-perturbative evaluation of the running coupling constant and a numerical check of Gribov's 
confinement scenarios for these two gauges. Our simulations are done in pure SU(2) lattice gauge theory at 
P — 2.2. We consider several lattice volumes in order to control finite- volume effects and extrapolate our results 
to infinite lattice volume. 



1. INTRODUCTION 

An essential step for understanding and ex- 
tracting physical information from gauge theories 
is the elimination of redundant gauge degrees of 
freedom. This is usually done by choosing a rep- 
resentative on each orbit of gauge-related fields 
(gauge fixing). In Ref. [111 Gribov showed that 
the Coulomb and Landau gauge-fixing conditions 
do not fix the gauge fields uniquely, i.e. there 
are many gauge-equivalent configurations satis- 
fying the Coulomb or Landau transversality con- 
dition. These Gribov copies do not affect pertur- 
bative calculations, but their elimination could 
play a crucial role for non-perturbative features 
of gauge theories, such as color confinement and 
hadronization. 

In order to get rid of the problem of spurious 
gauge copies, Gribov restricted jlj the physical 
configuration space to the region f2 of transverse 
configurations, for which the Faddeev-Popov op- 
erator is non-negative.^ This region is delimited 
by the first Gribov horizon, defined as the set of 
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J In the Coulomb case the transversality condition and the 
positiveness of the Faddeev-Popov operator are satisfied 
on each time slice. 



configurations for which the smallest, non-trivial 
eigenvalue of the Faddeev-Popov operator is zero. 
We now know that fl is not free of Gribov copies 
and that the physical configuration space has to 
be identified with the fundamental modular region 
I^H- Nevertheless, the region £1 is of interest in 
numerical simulations, since it is the space of con- 
figurations satisfying the usual lattice Coulomb or 
Landau gauge condition. 

The restriction of the path integral, which de- 
fines the partition function, to the region f2 im- 
plies a rigorous inequality [Q for the Fourier com- 
ponents of the gluon field. From this inequality, 
which is a consequence only of the positiveness of 
the Faddeev-Popov operator, it follows that the 
region Q is bounded by a certain ellipsoid E. This 
bound causes a strong suppression of the (un- 
renormalized) transverse gluon propagator D tr in 
the infrared limit (IJQ]. More precisely, it was 
proven H that, in the infinite- volume limit, D tr 
vanishes at zero momentum, although the rate of 
approach to 0, as a function of the momentum 
or of L, was not established. This is in marked 
contrast to the divergence in the infrared limit of 
the free massless propagator. 

Finally, because of entropy considerations , 
the Euclidean probability gets concentrated near 
the first Gribov horizon^ where the inverse of the 

2 This has been verified numerically in the minimal Landau 
gauge ||. 
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Faddeev-Popov matrix diverges. This causes an 
enhancement of the ghost propagator G(k) in the 
infrared limit Jl|||. This enhancement is a clear 
indication of a long-range effect in the theory that 
may result in color confinement. 

The confinement scenario is particularly simple 
in the minimal Coulomb gauge where the ghost 
propagator determines directly the Coulomb in- 
teraction 0J^]. In fact, in this case, confine- 
ment of color, i.e. the enhancement at long range 
of the color-Coulomb potential V(R), is due to 
the enhancement of G{k) at small momenta. At 
the same time, the disappearance of gluons from 
the physical spectrum is manifested by the sup- 
pression at k = of the propagator Z?jj (fc, fei) 
of 3-dimensionally transverse would-be physical 
gluons. Remarkably, V(R) is the instantaneous 
part of the 4-4 component of the gluon propa- 
gator D±i(x, t), and is a renormalization-group- 
invariant quantity [Q. Its Fourier transform V(k) 
may serve to define the running coupling constant 
of QCD by considering xog1;(\k\) = ^ = 

g%t D 44 (k), where x = 12N/(11N - 2N f ) has 
been calculated in fjj. Clearly, if the color- 
Coulomb potential V(R) is governed by a string 

~ -4 

tension at large distances, i.e. V(k) goes like 1/k 

at small momenta, then g^ik) ~ 1/k in the in- 
frared limit. 

Similarly, in Landau gauge, one can con- 
sider Q the running coupling constant 
gl(k) = g 2 [k 2 D(k)}[k 2 G(k)} 2 . This is also 
a renormalization-group-invariant quantity since 
(in Landau gauge) Z g Z^ 2 Z% = Z\ = 1. In 
this case we obtain g\{k) ~ k~ 2 if, for example, 
D tr (k) ~ const and G(k) ~ fc~ 4 in the infrared 
limit. On the contrary, if the gluon progator goes 
to in the infrared limit and the ghost propaga- 
tor blows up not faster than k~ A then g 2 L (k) has 
an infrared fixed point H], 

We test these theoretical predictions with data 
from a numerical study of SU{2) lattice gauge 
theory, without quarks, in the minimal Coulomb 
and Landau gauge at (3 — 2.2 . Simulations 
were done at different lattice volumes L 4 , in or- 
der to check for finite-size effects and, if possi- 
ble, to extrapolate to infinite lattice volume. De- 



tails of notation and numerical simulations are 
given, for the Coulomb case, in (h]] and will 
be presented, for the Landau case, in |11|. For 



these simulations we have used a cluster of AL- 
PHA work-stations (Coulomb and Landau data) 
and a PC cluster (Landau data) at the Dept. of 
Mathematical Physics of the University of Sao 
Paulo (DFMA/USP)f] In the Landau case we 
have also used a PC cluster^] at the Institute of 
Physics of the University of Sao Paulo, Sao Carlos 
(IFSC/USP). 
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Figure 1. Plot of the gluon propagators D tr (k) 
(lower curve) and D^k) (upper curve) as a func- 

-2 

tion of the square of the lattice momentum k 
for L — 28 (symbols * and O respectively) and 
L = 30 (symbols A and v respectively). Notice 
the logarithmic scale on the y axis. 



2. RESULTS 

Our data in the minimal Coulomb gauge |l(], 
clearly show (see Fig. |l|) that the equal- 
time transverse gluon propagator D tr (k, L) passes 
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Figure 2. Plot of the gluon propagator D tr (fc) 
(lower curve) and of the ghost propagator G(k) 
(upper curve) as a function of the square of the 
lattice momentum fc 2 for L — 22 (symbols * and 
O respectively) and L — 26 (symbols A and V 
respectively). Notice the logarithmic scale on the 
y axis. 

through a maximum and decreases as the momen- 
tum k approaches (for fixed L), and that V(k) 

-a -> 

is more singular than 1/fc at low fc, which in- 
deed corresponds to a long-range color-Coulomb 
potential. We have also obtained an excellent 
2-pole fit for D tT (k, L). Our fit indicates that the 
poles occur at complex m 2 = x(L) ± iy(L). In 
the infinite-volume limit we have m 2 — ± iy, 
for y = 0.375 ± 0.162 in lattice units, or y = 
0.330 ± 0.142 GeV 2 for the location of the gluon 

-.2 

poles in k . (It follows from the Nielsen identities 
|l3f that these poles are independent of the gauge 
parameters.) 

Similarly, in Landau gauge (see Fig. ||), the 
transverse gluon propagator D tr (k,L) is sup- 
pressed in the infrared limit, while the ghost prop- 
agator G(k) is more singular than 1/fc 2 . Also here 
a 2-pole fit can probably be used to fit the data 
for D tr (k,L). In this case, however, we need to 
improve the statistics in order to have better con- 
trol over the fits, and we need to simulate at larger 
lattice volumes (we went up to 26 4 ) to probe the 



infinite- volume limit. 

Finally, for the running coupling constant 
g 2 we consider jl0| the fitting formula k 2 = 
A 2 exp^g 2 )' 1 ] [{bg 2 y r + z{bg 2 ) a ]- 1 , which im- 
plicitly defines g 2 . Here r = 102/121, and A, 6, z 
and a are fitting parameters. For small g 2 , which 
corresponds to large momenta, this formula is 
dominated by the first term in the denomina- 
tor, whereas for large g 2 , i.e. small momenta, it is 
dominated by the second term in the denomina- 
tor. In particular, a governs the strength of the 
singularity of g 2 in the infrared limit; the case 
g 2 ~ 1/fc 2 corresponds to a = 1. 

In the Coulomb case we obtain [|l0| a = 1.9 ± 

0. 3. this corresponds to g 2 -, ~ l/|fc|, and V(k) ~ 
l/|fc| 3 at low momentum. In the Landau case we 
have a = 2.5 ± 0.3, which implies g 2 ~ 1/fc - 8 . 
In both cases we have fitted the data with a low 

->2 

momentum cut-off at k = 0.5. Of course, an ex- 
trapolation in (3 will be necessary to determine 
the strength of these singularities in the contin- 
uum limit. 
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